Universality in rotating strongly interacting gases 
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We analytically determine the properties of two interacting particles in a harmonic trap subject 
to a rotation or a uniform synthetic magnetic field, where the spherical symmetry of the relative 
Hamiltonian is preserved. Thermodynamic quantities such as the entropy and energy are calculated 
via the second order quantum cluster expansion. We find that in the strongly interacting regime the 
energy is universal, however the entropy changes as a function of the rotation or synthetic magnetic 
field strength. 

PACS numbers: 03.75.Hh, 03.75.Ss, 67.85-d 



Over the last few years ultracold degenerate gases have 
attracted much interest due to their controllability and 
stability. Advances in tight confining harmonic traps and 
the use of magnetic fields and Feshbach resonances in 
controlling atomic collisions have made it possible to ex- 
plore the BCS-BEC crossover [TH3]. Difficuhies with de- 
veloping a many-body theory for these systems in the 
strongly interacting regime using mean-field approxima- 
tions have motivated the study of few-body problems as 
a means to gain insight into the many-body problem. 
Few-body systems with contact interactions are exactly 
solvable or numericaly tractable [IHZj , particularly in the 
strongly interacting regime and have been experimentally 
studied in their own right [8] . The virial expansion of few- 
body physics can be used to calculate the thermodynam- 
ics of many-body systems [5HTT] and has been verified 
experimentally |12| . 

In this work we address the problem of unitary gases 
subject to a rotation or synthetic magnetic field by solv- 
ing the two-body problem and finding the virial expan- 
sion to second order. This enables us to show that en- 
tropy in the presence of a rotation or synthetic magnetic 
field is not universal, in contrast to the universal charac- 
ter of the total energy. 

A system subject to a rotation and one subject to a 
synthetic magnetic field have several similarities. In both 
systems angular momentum states and time-reversal 
symmetry are broken. Furthermore, both problems can 
be described by gauge-dependent Hamiltonians, making 
it convenient to consider the systems together and to 
draw comparisons between the two. In ultracold trapped 
gases the dominant contribution to the low energy be- 
havior is from the two-particle s-wave interactions. 

To begin the analysis the rotating system is considered 
first. Specifically, the motion of two particles of mass m 
in a harmonic trap potential Vtrap(r) subject to a rotation 
Q, and a contact interaction potential Vint(i"i — are 
described by the Hamiltonian 



H 



P? 



Pi 



Vtrap(ri) + Vtrap(r2) + 



2m 2m 

rj • ri X pi + rj • ra X p2 + l^„t(ri - ra), (1) 

where rj and pi are the positions and momenta of each 
particle. Equation ([T]) can be decoupled in center of mass 



and relative coordinates, yielding 
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ytrap(R)-f^-RxP (2) 
Vtrap(r)-nTXp + l^nt(r), (3) 



where R = (ri + r2)/\/2 and r — (ri — r2)/V2 arc the 
center of mass and relative coordinates. We consider the 
case where the rotation is about the z-axis with frequency 
ilz so that rJ = (0,0, r^z). The harmonic trapping po- 
tential is chosen to be axially symmetric with transverse 
and axial frequencies uj_i and w^, respectively. Using the 
axial trap length d ~ y/h/{mujz) and energy hujz , Eq. ^ 
can be written in the dimensionless form 
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where rj = uj^/uz and ^ — —Slzl^z parameterize the as- 
pect ratio of the trap and the applied rotation. Equation 
(|4]) is the Hamiltonian for a shifted anisotropic harmonic 
oscillator with eigenstates 
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and eigenenergics, in units of fvjjz, 

Enmk ^ {271 +\m\ + l)7j + m^+{k + 1/2), (8) 

where l}^\p) and Hk{z) are, respectively, Laguerre and 
Hermite polynomials. The center-of-mass component is 
therefore solved exactly and the effects of the interparti- 
cle interaction are described entirely within the relative 
component. 

In the presence of a synthetic magnetic field the 
Hamiltonian for the center-of-mass of two anisotropically 
trapped particles is 
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Krap (R) , 
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where A is the synthetic magnetic vector potential. 
While the properties of a magnetic system are formally 
independent of choice of gauge, the problem is not ana- 
lytically tractable in the Landau gauge, but is tractable 
in the symmetric gauge. For a uniform synthetic mag- 
netic field, B, in the z-direction wc therefore set A = 
^{—y,x,0). It can then be shown that the eigenstates 
and eigenenergies for Eq. ^ are given exactly by Eqs. ([s]) 
and Q with redfined parameters 



and r]c 
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FIG. 1. Eigenenergies of the relative Hamiltonian for (a) 77 
1 and (b) 77 = \/ll/9 
scattering length. 



1.10, as a function of the inverse 



where ujc — \qB/M\ is the cyclotron frequency with 
M = 2m. For the single-particle Hamiltonian, the eigen- 
states and eigenenergies are expressed in terms of f 1 and 
771 using M — m and for the relative Hamiltonian the 
eigenstates and eigenenergies are expressed in terms of 
^roi and 77i-ci using M = m/2. In contrast, for the rotating 
case the relative, centre-of-mass and single particles are 
parameterised by the same 77 and ^ as defined in Eq. Q . 

Having solved the center of mass Hamiltonian in both 
the rotating and synthetic magnetic field cases, the rel- 
ative Hamiltonian needs to be solved exactly. The low- 
energy regularized s-wave contact interaction is |13) 



^nt(r) 
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where a is the scattering length. In the single-particle 
basis the relative wavefunction can be written as 
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where the m ^ states are omitted because they do not 
contribute, due to the centrifugal barrier in cylindrical co- 
ordinates. As such, the projection of angular momentum 
along the z-axis due to the rotation or synthetic mag- 
netic field is contained entirely within the center-of-mass 
energy. Following [141 we find that the energy spectrum 
of the relative motion can be determined from 
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This integral is not analytic in general and is formally 
valid only for a; > but can be extended to all energies 
by the recurrence relation 
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In Fig. [T] we plot the energy eigenspectrum Eq. ( 13 ) 
as a function of inverse scattering length d/a for two val- 
ues of 77. For the case 77 = 1, Fig. [ija) shows states 
evenly spaced by 2hiOz with a smooth transition across 
the Feshbach resonance and a single bound state as the 
ground-state in the repulsive regime. Unlike in the spher- 
ical basis, the s-wave interaction allows the orbital angu- 
lar momentum states to be present at each energy level 
but they are degenerate. It is possible to take a linear 
combination of these degenerate wave functions to obtain 
the usual spherically symmetric, zero orbital angular mo- 
mentum wave function. If 77 is rational, by the properties 
of Q some of these degeneracies are restored. For the 
special case when 77, or I/77, is an integer then all de- 
generacies are restored [T3]. The most extreme case is 
when 77 is irrational, for which there can be no degen- 
eracies with higher states. If the spherical symmetry is 
broken so that 77 7^ 1 then the degeneracies are lifted, 
as shown in Fig. [ijb). Angular momentum is no longer 
conserved and a linear combination of wave functions will 
not reduce the number of accessible states. As such, one 
could argue that as soon as the spherical symmetry of 
the problem is broken the properties of the gas become 
independent of the relative Hamiltonian and hence the 
interactions. Of course, this cannot be correct and to 
overcome this dichotomy one must consider an interac- 
tion potential that is regularized in accordance with the 
symmetries of the system |16) . and including contribu- 
tions from higher order (non-s-wave scattering) processes 

[IIHID]. 

In the following analysis, we assume spherical symme- 
try of the relative Hamiltonian. In the rotation case, 
r] — 1 means that the trap is isotropic and the rela- 
tive eigenenergies are unaffected by the rotation. In the 
synthetic magnetic field case, 771.0I = 1 is equivalent to 
w_L = [wf — (gS/rTz)^]^/^, fromEq. (10). This means that 
for a given magnetic field there is a particular trap shape 
that preserves the spherical symmetry of -ffrci- 

Few-body physics has importance beyond small scale 
systems into the thermodynamics of a many-particle gas. 
We can achieve this through a quantum cluster expansion 
of the grand thermodynamic potential fl = —ksTlnZ in 
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terms of the fugacity z = exp{fi/kBT): 
n = -keTQi {z + b2Z^ + ■ 

where 

62 = iQ2 - Ql/2)/Qi 



(16) 
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is the second virial coefficient and the iV-particle parti- 
tion function Qn = Tr[exp(— ^^jv/ZcsT)] is determined 
from solving the A^-body problem [Tlj . 

In order to calculate the thermodynamics it is more 
convenient to consider the difference between the inter- 
acting and non-interacting systems and define 



Ab2 = {b2-b'^^) = iQ2- 



(18) 



where the superscript '0' denotes non-interacting quan- 
tities. The thermodynamics of a non-interacting two- 
component Fermi gas in an anisotropic trap with a con- 
stant magnetic field or a rotation are determined from 
the grand potential fi'^^'. Using the energy spectrum Q 
to determine the density of states it can be shown that 
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Equations p9| and (20) reduce to the more familiar 



forms in the limits of no rotation or synthetic magnetic 
field {£_ 0) and isotropic trapping (77 1) [10]. The 
integral in Eq. ( 19 ) can be expanded in powers of the fu- 



gacity z to obtain the non-interacting virial coefficients 



In the thermodynamic limit, A62 is independent of 
temperature. In the following calculations it is useful to 
introduce the small parameter a) = fiw^/fesT. Expand- 
ing the virial coefficients and cluster partition function 
in the small parameter ui allows us to determine their 
universality. 

Since the two-body problem may be separated into rel- 
ative and center-of-mass coordinates, we can sum over 
the center-of-mass component Qcm and the relative com- 



ponent independently so that Eq. ( 18 ) becomes 
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In the case of a rotating gas both Qi and Qcm are 
determined from the single particle energy spectrum (|8]). 
Including a factor of 2 to account for the spin states 
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and Qcm = Qi/2. In the high temperature limit ((D — 0) 
the leading order behavior of Eq. ( 22 1 is exactly Eq. ( 20 1 . 



The magnetic field case is more complicated. As in the 
rotating case the same energy spectrum is used. However, 
Qi is obtained by exchanging 77 and ^ in Eq. (22) with 
771 and ^1. Similarly, Qcm is obtained by exchanging 77 
and with 77cm and ^cm, and omitting the spin-counting 
factor of 2. 



To perform the remaining sums in Eq. ( 21 ) we need the 
eigenenergies of -ffrci. Specifically, the case of 77 = 771.01 = 
1, i.e. the relative Hamiltonian is isotropic, is considered. 
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The non-interacting (a = 0) spectrum is 
3/2 and the spectrum in the unitary regime {a ^ 00) is 
i^roi = 2n -I- 1/2 4. For attractive interactions all states 
are included, but for repulsive interactions the 71 = 
bound state is omitted. 

For a rotating trapped gas in the high temperature 
limit with a large number of particles 
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and for a trapped gas in a synthetic magnetic field 
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We are now able to calculate thermodynamic quantities 
like the total energy and entropy of the gas from the 
grand potential Eq. (16) 



To determine the energy and entropy the fugacity z 
must be calculated first from N{z) — —dfl/d^, which 
is quadratic in z. In the rotating case there is always 
a single positive root of N{z) and so the fugacity is al- 
ways physical. In contrast, in the synthetic magnetic field 
case the roots are non-trivial. Despite the temperature 
regime relative to the trap ground-state energy being set 
by uj, the temperature decreases as ^1 icreases, where ^1 
parameterizes the synthetic magnetic field. This is due 
to the restriction of isotropy in the relative coordinate, 
i.e. rj^ei = 1 so that the transverse trapping must change 
as the magnetic field changes. For the attractive case, 
in all regimes there is a pole at ^1 — 4/\/3T « 0.718. 
For magnetic field strengths below this there is always 
one non-negative solution, which is physical. Above this 
value, we see that for larger uj more critical points appear 

at = 4/^/28 w 0.756 and Ci = 4/%/T9 « 0.918. Specif- 

,1/3 



ically, this occurs when w > [(820 + 32V65T)/(9iV)] ' 
Between these two critical points there is always one 
physical solution, but outside these values both solutions 
are complex. As a result of this behavior, we consider 
that the virial expansion for the synthetic magnetic field 
case in the attractive regime is valid for ^1 < 0.75. This 
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FIG. 2. Energy per particle (a), (b) and entropy per particle 
(c), (d) of a rotating atomic gas and a uniform synthetic mag- 
netic field, respectively, for attractive interactions. Curves are 
plotted for 5 = 5i = 0.0 (solid), 0.5 (dashed), 0.75 (dotted) 
for both systems and additionally for ^ = 0.95 (dot-dashed) 
for the rotating system (a,c). 

problem is not present for the rotating gas where the 
parametrization merely restricts ^ < 1, which is not an 
issue because we already expect that when \^z\ > 
the gas becomes effectively untrapped. The same cal- 
culations can be performed for repulsive interactions to 
determine where the virial expansion is valid. 

For the case of attractive interactions, we plot in Fig. [2] 
the energy (a,b) and entropy (c,d) per particle of a gas in 
a harmonic trap, subject to a rotation (a,c) or a synthetic 
magnetic field (b,d) in the case rj — rj^-d = 1. Increasing 



the rotation or the synthetic magnetic field at a given 
temperature does not significantly affect the total energy 
of the system. In both cases the rotation and synthetic 
magnetic field can be viewed as a reparameterization of 
the transverse trapping frequencies, but at unitarity the 
interactions still dominate the energy of the system. In 
contrast, for both systems the entropy increases for larger 
^ (or ^i), as more states become accessible to the centre- 
of-mass because the rotation or synthetic magnetic field 
couple higher angular momentum states to lower ener- 
gies. In the case of repulsive interactions, although the 
details are different, the same calculations can be per- 
formed and the results are qualitatively the same; the 
total energy of the gas is universal, but the entropy in- 
creases as rotation frequency or synthetic magnetic field 
strength is increased. 

In conclusion, we have parameterized the problem of 
two ultracold atoms in a harmonic trap subject to a ro- 
tation or a synthetic magnetic field to retain spherical 
symmetry in the relative Hamiltonian. When the spheri- 
cal symmetry is broken, even by a small perturbation, the 
s-wave contact interaction allows many more distinct rel- 
ative energy states. These energy states appear because 
the spherical nature of the regularized interactions is in- 
compatible with the cylindrical symmetry of the trap. 

In the special case of the relative Hamiltonian being 
isotropic in the unitary limit the total energy in the ther- 
modynamic limit is universal and the entropy is not. This 
is due to the rotation or synthetic magnetic field cou- 
pling higher angular momentum states to lower energy 
levels and increasing the number of states available to 
the centre-of-mass of the system in a given energy range. 
However, in the unitary limit the total energy is domi- 
nated by the relative energy, which is determined by the 
interaction energy, and hence is largely independent of 
the rotation or synthetic magnetic field. 
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